A simple design of an artificial electromagnetic black hole 
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We study the properties of an artificial electromagnetic black hole for transverse magnetic modes 
rigorously. A multi-layered structure of such a black hole is then proposed as a reduced variety 
for easy experimental realizations. An actual design of composite materials based on the effective 
medium theory is given eventually with only five kinds of real isotropic materials. The finite element 
method confirms the functionality of such a simple design. 



I. INTRODUCTION 

Transformation opticsii^ is a useful tool to design many 
novel wave manipulation devices^^— . Its general form, 
the "general relativity in electrical engineering"—, was 
proposed to mimic cosmic optical properties^^*^. A fur- 
ther method was also suggested^^, which can transmute 
the anisotropic parameters from the general form into 
isotropic ones. An isotropic "optical black hole" was 
given as well based on such a transmuting methodic. 
Another approach to a broadband absorber was also pro- 
posed from Hamiltonian optics^^. Such an absorber can 
be termed as an effective 'optical black hole"—, which 
was later implemented by using nonrcsonant metama- 
tcrial units in microwave frequencies for transverse elec- 
tric(TE) modesii. The experiment demonstrated the im- 
portance of metamaterials in implementing the related 
transformation media devices once again. However, the 
structure is a bit complicated and will be very challenging 
to adapt for higher frequencies. 

In this paper, we will consider the transverse mag- 
netic(TM) modes of such an artificial electromag- 
netic (EM) black holeii and give a simple design with 
several kinds of real materials. In the following, we begin 
with the rigorous studies of TM modes of such an EM 
black hole in section 2. Then we propose a simple design 
based on composite materials in section 3. Finally, we 
give conclusions in section 4. 



II. THE RIGOROUS CALCULATIONS OF AN 
ARTIFICIAL BLACK HOLE FOR TM MODES 



We start from the same permittivity profile of the ar- 
tificial EM black hole as in Ref J^, 
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where Ri = Riy -p- with the geometry parameters de- 
scribed in FiglUa). 

As it is mentioned in Ref.— , the TE and TM polar- 
izations decouple and can be solved independently. The 



TE EM black hole has been studied rigorously^^ and im- 
plemented with nonrcsonant mctamaterial units^i. With 
similar procedures, we will focus on the TM modes with 
the magnetic field along the z-direction and propose a 
much simpler design in the following sections. 

In the region i?i < r < R2, the general solutions of the 
magnetic field in z-dircction can be written asi^i^^. 
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where n is the angular momentum number, u is the an- 
gular frequency, and HY^{r) satisfies. 



ir^M + {klRl n')Hr{r) = 0, (3) 



dr^ dr 
whose solutions can be expressed as 



HTir)=A„[ — 
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with the wave vector of light in vacuum fcp = ^ . 

In the region r > R2, we suppose that an Gaus- 
sian beam is incident with its magnetic field in z- 
directioi*i2^, 
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where Jn{kor) are Bessel functions of order n, and P„ is 
determined hy^^^, 
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where A is the wavelength of the Gaussian beam, Wq is 
half of the beam waist, (a;o,?/o) is beam center in Carte- 
sian coordinates, and a = sin~^(-^) is the angle between 
the wave vector and i-axis. The scattering waves are 
assumed to be, 
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where H„ {kor) are the first kind Hankel functions with 
the angular momentum number n. Hence the total mag- 
netic field 7?°"* (r, 6) is the summation of the incident 
field Hl'"={r,9) and the scattering field iJf™(r,e'), 

OO 
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X exp(m0 — itjt). 

In the region r < Ri, the absorbing core is isotropic 
with permittivity e'+ie" , the magnetic field in z-direction 
can be written as, 

OO 
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(9) 

From the continuous conditions on the boundaries r = 
Ri and r = R2, i.e., the continuities of the magnetic field 
in z-dircction and its normal derivative, the coefficients 
Qn, Bn, and C„ can be uniquely determined in terms 
of the incident coefficients P„ . 

In this paper, we will set i?i = 60mm, R2 = 120mm, 
e' = 4.0, e" = 0.33^, and A = 30?72m for instance. The 
waist of the incident Gaussian beam is 60mm and the 
beam center is at (—120mm, 60mm). We use a cut-off an- 
gular momentum number rimax = 100 during the calcula- 
tions. The magnetic field intensity pattern |i?zp is plot- 
ted in FigHJb) , which shows that the beam is bent enor- 
mously toward and absorbed by the inner core. Thereby 
the above absorbing system can be termed as an "EM 
black hole"i^. 

The absorption cross section per unit length can be 
written as2^, 
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with Dn = Qn/Pn- Substituting the above parameters, 
the absorption cross section of the present artificial black 
hole is 0.2231, which is close to its geometry cross section 
2i?2^. That means this artificial black hole can sever as 
a nearly perfect absorber as in Ref J^. 

III. TWO-STEP APPROACH TO A SIMPLER 
DESIGN OF AN EM BLACK HOLE, A LAYERED 
STRUCTURE AND AN ACTUAL DESIGN 

The above artificial black hole can be designed by us- 
ing multi-layered cylindrical structuroi^. We break up 
the inhomogeneous region Ri < r < R2 into N concen- 
tric shells of isotropic dielectrics with equal thickness. In 
FiglUa), we plot the relationship between the layer num- 
ber N in Ri < r < R2 and the absorption cross section, 
which is calculated by generalized Mie theory (GMT) for 
multi-layered structure and Eq. (|10p . The results show 
that a twelve-layer structure is good enough to imple- 
ment such an EM black hole. As a concrete example, 




FIG. 1: (a). The schematic plot of an incident Gaussian beam 
interacting with an artificial EM black hole. R2 is the radius 
of the "event horizon" described in Ref.— , _Ri is the radius 
of the inner absorbing core. (b). The magnetic field intensity 
pattern for an incident Gaussian beam interacting with the 
EM black hole, which is calculated based on rigorous solutions 
Eqs. Q-®. 



we will use 12 layers of isotropic materials whose rela- 
tive permittivities are shown in Figl^l^b). The outmost 
layer is set to be air while the permittivity of the inner 
absorbing core is Ecore — 4: + 0.33»2i. Figl^^c) shows the 
magnetic field intensity pattern near the present layered 
structure black hole based on GMT. The related absorp- 
tion cross section is about 0.2222, which means that the 
layered structure works as well as the original one. 

However, it would be difficult to find the above twelve 
layers of isotropic materials one by one. Based on a two- 
step approach^^, we will implement this layered EM black 
hole with only five kinds of real isotropic materials. We 
will use the following real materials, air, aluminum (Al) 
metal rods, polyethylene (PE), polymethyl methacry- 
latc (PMMA) plexiglass, and Polyvinylidene fiuoride 
(PVDF). Their related permittivities in microwave fre- 
quencies (near 10 GHz) are about epE = 2.3, epMMA = 
2.6, and epvDF = T + 0.7»2i (Note that fair = 1 and 
tAi — —00^). As the equal thickness of each layer is 5mm, 
we divide each layer into cells each measuring 5mm in the 
middle so that each cell has "fanlike" shaped. A cylin- 
der with one kind of materials is placed in the center of 
each cell and embedded in a background material. Be- 
cause the thickness of each cell/layer is very small when 
compared with the wavelength, we can approximate the 
above twelve layers of isotropic materials with compos- 
ite materials of the above five kinds of materials based 
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FIG. 2: (a). The absorption cross sections varies with dif- 
ferent layer numbers in a layered EM black hole. (b). The 
relative permittivity for each layer (red lines) in a twelve-layer 
EM black hole. The black curve is e = {R2/r)^ for compari- 
son, (c). The magnetic field intensity pattern near the layered 
EM black hole, calculated with GMT. 



on the effective medium theory (EMT)^"— . Each layer 
of isotropic materials can be approximated by a circular 
array of cylinders embedded in a background material, 
see for details in Figl3Ja). From the EMT, the effective 
permittivity e^ff of composite materials with cylinders 
embedded in a background material with square lattice 
satisfies, 



ei - e 



eeff + e ei+e 



fs J 



(11) 



where e is the permittivity of the background material 
while ei is the permittivity of cylinders with square lat- 
tice, and fs is the filling ratio of the cylinders. Suppose 
the lattice constant is a = Bmm, fs = with the 

radii of cylinders re- 

As the permittivities of the outer nine layers are in 
the range of 1 and 2.3, we can use air-PE composite ma- 
terials to approximate each layer. We shall embed air 
hole cylinders in PE to approximate the permittivities 
from 1.28 to 2.12 (seven layers). The radii of cylinders in 
each layer can be obtained from Eq. ([TT|) (from 2 Aram 
to 0.9mm, see also in Figl2l^a), the PE is denoted with 
green color while the air hole cylinders are denoted with 
blue color). As permittivities of the inner three layers in 
i?i < r < i?2 is 3.66, 3.15 and 2.6, we can use Al-PMMA 
composite materials to approximate two of them. As the 
effective permeability of such composite materials is not 



unity2^, we shall use the square of the effective refractive 
index (but not the effective permittivity) to approximate 
the above two permittivities (3.66 and 3.15), i.e., 



T-lff = t^<^eff =e{l + fs). 



(12) 



Likewise, we can obtain the radii of Al cylinders (1.8mm 
and 1.3mm from Eq. ((T2|) . see also in Figl3l[a), the 
PMMA is denoted with orange color while the Al cylin- 
ders are denoted with dark blue color). The inner absorb- 
ing core can be approximated by using air hole cylinders 
embedded in PVDF. The inner core is also divided into 
twelve layers with the same procedure as that applied 
to the region i?i < r < R2. The radii of the air hole 
cylinders are 1.7mm from Eq. ([TT|) (see also in Figl3l[a), 
the PVDF is denoted with brown color while the air hole 
cylinders are denoted with blue color). 
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FIG. 3: (a). An actual design of the EM black hole with com- 
posite materials. The air, Al, PE, PMMA, and PVDF are 
denoted by blue, dark blue, green, orange, and brown colors 
respectively. The radii of cylinders in each layer (vc) are also 
given correspondingly, (b). The magnetic field intensity pat- 
tern near the black hole with composite materials described 
in (a), (c). The same to (b) but replacing the inner absorb- 
ing core with PVDF directly. The simulation results are from 
COMSOL Multiphysics finite element-based electromagnetics 
solver and the incident beam is a fundamental-mode Gaussian 
beam for approximations to Eqs. ((5)1 and (j6]). 



Figure |31[b) shows the magnetic field intensity pattern 
near the black hole designed above in Fig|31[a). We use 
COMSOL Multiphysics finite element-based electromag- 
netics solver to perform the simulation, which shows that 
the black hole with the above composite materials can ap- 
proximate both the layered structure black hole and the 
original EM black hole very well. In addition, we found in 
FigEIc) that if we replace the inner core by using PVDF 
solely (without embedding the air holes), the device still 
functions. Hence these two designs both serve for easy 
experimental fabrication of EM black holes. 
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IV. CONCLUSION 

To summarize, we have rigorously studied the proper- 
ties of an EM black hole for TM modes. With a two-step 
approach, we proposed actual designs of such a black hole 
by using composite materials with only five kinds of real 
isotropic materials, enabling easy fabrications with nowa- 
days technologies. Our design circumvents retrieving spe- 
cific constitute parameters from resonant structures, the 
device is expected to function in a broad bandwidth of 
frequencies. In particular, due to the simplicities of such 
designs, it would be feasible to adapt the same method 



for higher frequencies, such as THz, infrared, or even op- 
tical frequencies^!^. 
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